Oblate-prolate shape coexistence at finite angular momentum 
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We investigate shape coexistence in a rotating nucleus. We concentrate on the interesting case 
of ^■^Kr which exhibits an interesting interplay between prolate and oblate states as a function of 
angular momentum. The calculation uses the local harmonic version of the method of self-consistent 
adiabatic large-amplitude collective motion. We find that the collective behaviour of the system 
changes with angular momentum and we focus on the role of non-axial shapes. 

PACS numbers: 21.60.-n, 21.60.Jz 
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Only a small number of quantum mechanical many- 
body systems can be solved exactly through analytical 
or numerical means. Many approximate schemes exist 
to find reasonable answers, but these may not always be 
easy to interpret in physical terms. Another slant on this 
problem is to try and describe some of the dynamics of 
such systems in terms of a limited set of degrees of free- 
dom. These should of course be chosen through some 
method appropriate for the problem at hand. Many ap- 
proaches are available, in areas ranging from field theo- 
ries to atomic physics (see, e.g., the reviews in Ref. [l|). 
These are typically based on the concept of "relevant de- 
grees of freedom" , or on the introduction of collective mo- 
tion and collective paths - which are two ways to express 
rather similar principles! Depending on the energy scales 
involved and the physical situation being described, the 
remaining degrees of freedom are either frozen out, or 
treated as a heat-bath for the motion of the relevant de- 
grees of freedom. 

The most common description of nuclear collective dy- 
namics is influenced by the success of the liquid drop 
model _2J. The description is based on the competi- 
tion between quadrupole shape fluctuations of the nu- 
cleus, typically modelled microscopically by a long-range 
quadrupole-quadrupole force, and a BCS-like pair con- 
densate, modelled by a short-range pairing force. The 
model of the nucleus does not have to be this naive, but 
even with better models, it is often assumed that the low- 
energy dynamics is well described by quadrupole shape 
fluctuations. We would like to question this assumption, 
which leads to a search for the best choice of collective 
coordinates (the standard terminology for the relevant 
degrees of freedom in this field). There are quite a few 
partial answers, see the review [3| for a discussion of some 
of these. The goal is to find a method that determines a 
collective path self-consistently, based only on knowledge 
of the Hamiltonian governing the system. 

Clearly, in nuclear physics one does not know the 
Hamiltonian. However, for the collective properties of 
medium mass to heavy nuclei, models that contain the 
two key parts of the nuclear force, a short-range pairing 
force, and a long range multipole-multipole force (usu- 
ally approximated by a quadrupole-quadrupole one) are 



known to be able to capture the essential part of the 
physics, see e.g. the textbook Ref. |j|. 

With modern experimental techniques, we can create 
nuclei at considerable angular momentum, and study the 
behaviour as a function of angular momentum, which 
provides us with an additional parameter we are almost 
free to choose. The resulting experimental data are of- 
ten analysed in terms of simple collective models, usually 
expressed in terms of quadrupole shapes. The physics 
of collective motion is then often described in terms of 
shape transitions and shape coexistence, the mixing be- 
tween various shapes. One important example is the 
recent large interest in shape coexistence at low angu- 
lar momentum in nuclei in the A ~ 70 and N k, Z 
mass-region |a, Ifl LD, Isl Isl- In this region one finds a 
large diversity of shapes, and rapid changes in shape 
with particle number and angular momentum. This large 
variety is caused by shell gaps corresponding to spheri- 
cal {N,Z = 38), prolate (N,Z = 34,38,40) and oblate 
(N , Z — 34, 36) shapes that exist in this mass-region. 
The detailed understanding of how states with differ- 
ent shapes mix is, among others, important for the as- 
trophysical nucleosynthesis process that passes through 
these proton-rich nuclei. 

The quality of the recent experimental data, fuelled by 
advances within the field of 7-ray spectroscopy, has led 
to theoretical effort focused on this mass-region. The nu- 
cleus ^^Kr, the subject of this letter, shows oblate-prolate 
shape coexistence and/or shape transitions, the nature of 
which seems to depend strongly on angular momentum. 
There are other calculations for this nucleus: Ref. ^] have 
used the constrained mean-field method to study the po- 
tential energy surface. The approach used in Ref. [2| goes 
beyond the mean field approximation but does not an- 
swer the question which degrees of freedom are important 
for the collective path. Oblate-prolate shape coexistence 
has also recently been studied in other mass regions, see 
e.g. Ref. [13. 

Our formalism, as set out in detail in [3,lll|, is based 
on time-dependent mean field theory, and the fact that 
a classical dynamics can be associated with it. The is- 
sue of selecting collective coordinates, and determining 
their coupling to other degrees of freedom, thus becomes 



an exercise in classical mechanics. If we assume slow 
motion the mean-field energy, which is also the classi- 
cal Hamiltonian, can be expanded to second order. This 
corresponds to a parametrisation of the one-body den- 
sity matrix in terms of a set of canonical coordinates, 
f", and conjugate momenta, tTq,. The potential V{^) and 
the mass matrix B°'^{£) are the coefficients of the ex- 
pansion of the classical Hamiltonian in powers of the mo- 
menta TT at zeroth and second order, respectively. Within 
this adiabatic Hamiltonian we search for collective (and 
non-collective) coordinates q^ and conjugate momenta 
Pp. These are assumed to be obtained by an invertible 
point transformation of the original coordinates ^" and 
momenta tTq , preserving the quadratic truncation of the 
momentum dependence of the Hamiltonian, by 

g^^riO, Pp=3>a (m, a =1, ••.,«) (1) 

where we use a standard notation for the derivatives. The 
adiabatic Hamiltonian is then transformed into 
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n{q,p) = v{q) + -B^^y^p, + o(/) 



(2) 



in the new coordinates. The new coordinates q^ are now 
to be divided into three categories: the collective coordi- 
nate, the zero-mode coordinates, which describe motions 
that do not change the energies and finally the remaining 
non-collective coordinates. 

In the local harmonic approximation (LHA) the col- 
lective coordinate is determined by means of the solution 
to a set of self-consistent equations. These are the force 
equations and the local RPA equation 



H.. 



Kf^a+t^Jx,a+ ^ Mt^^t.q, (3) 



^^B-^^fji = {hnff^c. 



(4) 



The parameters u and /i are Lagrange multipliers that 
implement the condition of fixed angular momentum 
along the x axis, and fixed particle number, respectively. 
In nuclear physics these are usually called generalised 
cranking parameters. The parameter A is an additional 
Lagrange multiplier for the collective mode, forcing the 
system to stabilise at a point away from equilibrium. The 
covariant derivative Va/s is defined in the usual way [13 . 
The collective path is found by solving Eqs. Q and Q 
self-consistently, i.e., we look for a path consisting of a 
series of points where the lowest non-spurious eigenvector 
of the local RPA equations also fulfils the force condition. 
In the minimum of the potential the spurious solutions 
decouples from the other collective and non-collective so- 
lutions. In this paper we have chosen to ignore the effects 
of the spurious admixtures to the RPA wave- functions, 
but test calculations has showed these to be small. To 
limit the computational effort we use the method pre- 
sented in Ref. [llj, [l^ to reduce the size of the RPA 
matrix. There it was shown that the RPA equation can 




FIG. 1: Large amplitude collective motion in ^^Kr at 7 = 0. 
Q is the collective coordinate, a) Energy along the collective 
path, b) The square of the lowest RPA frequencies, c) The 
deformation £2 and the triaxiality 7. d) The dimensionless 
pairing operators (Pt)- 



be solved with good accuracy by assuming that the RPA 
eigenvectors can be described as a linear combination of 
a small number of state-dependent one-body operators. 
The same basis set turned out to work well also at finite 
angular momentum. 

We apply the LHA to the constrained pair- 
ing-|-quadrupole Hamiltonian as described in [l3| with a 
constraint on particle numbers and angular momentum 



H' ^ h^-Y^GrPlPr-l Y. QmQm-^J. 



Z^ 



M=-2 



UrNr, 

(5) 

where ho is the spherical Nilsson Hamiltonian jlj]. Qm 
and Pr are the dimensionless quadrupole and pairing op- 
erators 13]. This Hamiltonian is treated in the Hartree- 
Bogoliubov approximation, within a model space consist- 
ing of two major shells. We follow Ref. [ijl and multiply 
all quadrupole matrix elements with a suppression factor. 
We start by examining the non-rotating states in ^^Kr 
where a prolate-oblate shape coexistence is established 
experimentally [5|. We find a collective path going from 
the oblate minimum over a spherical energy maximum 
into a prolate secondary minimum and continuing to- 
wards lager deformation, see Fig. ^ This is in contrast 
to the result recently found for ^^Se Q where the non- 
rotating prolate and oblate minima are connected via 
non-axial states. At a value of the collective coordinate 
Q ^ 7.5 there is an avoided crossing between the lowest 
RPA mode, which we are following, and two higher lying 
modes that are of pairing vibrational character as can be 
seen in Fig. ^. At this point our algorithm is no longer 
able to find a stable solution. This indicates that more 
than one collective coordinate should be used, which at 
the moment is a very difficult calculation. At large oblate 
deformation we see a collapse in neutron pairing after an 
avoided crossing between the /3-vibration we are follow- 
ing and a pairing vibration. As explained in detail in 





FIG. 2; Large amplitude collective motion in ^^Kt at / 
See Fig. for an explanation of the various figures. 



Ref. [ll| , such points are like the origin of radial coordi- 
nates for the collective motion, and it is not correct to 
continue the calculation. Clearly this does not preclude 
other collective coordinate becoming important here. 

At finite angular momentum the collective path will 
no longer go through the spherical state, since the rota- 
tional energy diverges for such a state. Instead the oblate 
and prolate minima are connected by a path consisting 
of non-axial states. Due to problems with narrow level 
crossings we have to start in both the prolate and oblate 
minima to find the complete collective path. In Fig. [21 
we see that if we start in the (almost) oblate minimum 
and follow the collective path towards larger oblate defor- 
mation the situation is very similar to the non-rotating 
case discussed above: We first see a decrease in the the 
quadrupole deformation, but after an avoided crossing 
with a pairing vibration the collective path changes its 
nature and ends in a neutron pair-field collapse. When we 
follow the collective path in the other direction, towards 
smaller oblate deformation the system goes through an 
avoided crossing with a 7-vibrational mode. The path 
then turns into the triaxial plane and the 62 deformation 
is almost constant, but with 7 decreasing quickly from 
60° to 0°. At this point the collective path goes through 
another avoided crossing this time with a /3-vibration. 
After the crossing the collective path follows an almost 
prolate shape with increasing deformation. The number 
of avoided crossing suggests that our description, even 
though it captures much of the physics, is not fully quan- 
titative. 

In Fig. 121 at angular momentum I — A we see a very 
different character of the collective path. The collective 
path is limited by the collapse of both the proton and 
neutron pair-field. The points of zero neutron or pro- 
ton pairing has an excitation energy of less the 150 keV. 
When starting the calculation in the prolate minimum 
we only find the collective path in a limited region due 
to avoided crossings. This suggests that this part of the 
collective dynamics is irrelevant. A related calculation 
at zero gap, which hardly changes the minimum energy. 



FIG. 3: Large amplitude collective motion in ^^Kr for I — A. 
Q is the collective coordinate, a) Energy along the collective 
path, b) The dimensionless pairing operators (Pr). 



but has different fluctuations, suggests that there is no 
collective dynamics around the 4+ state calculated here, 
and thus no shape mixing or shape coexistence. 

We want to solve the collective Hamiltonian along the 
collective path. After having made a semi-classical ap- 
proximation, which leads to a classical Hamiltonian, we 
need to remember that we are studying a quantum sys- 
tem. The standard technique to deal with this is to treat 
the classical Hamiltonian as a quantum one, and to cal- 
culate the eigenfunctions and energies. Details on how 
we do this can be found in Ref. |ll|. We have solved 
the collective Hamiltonian for the / = 0, 2 and 4 cases 
in ^^Kr. We also need to calculate the proton and neu- 
tron pairing-rotational masses. They becomes zero at 
the point where (P) — > 0. In Fig. 21 we see that the 
ground state wave-function is concentrated in the oblate 
minimum, but that it has a substantial spread along the 
collective path and is skewed towards the spherical state 
due to the collapse in the neutron pair-field. The first ex- 
cited / = state has its major component in the prolate 
minimum with a small component in the oblate mini- 
mum. The prolate peak in the wave-function for the first 
excited state is broader and more symmetric than for 
the ground state. The third / = state is approximate 
spherical but has substantial oblate as well as prolate 
components. The prolate state lies at a very low excita- 
tion energy of only 0.37 MeV, while the spherical state 
is somewhat higher in energy. For 1 = 2 the situation is 
similar to the case / = 0. However, the mixing of the pro- 
late and oblate states is substantially stronger. The long 
tail of the oblate peak stretches along the collective path 
into the prolate minimum, and that there is a secondary 
oblate peak in the prolate state. Due to the wide peaks 
of the collective wave function the expectation value of 
deformation in the collective states are substantially dif- 
ferent from those of the minimum potential energy states, 
as can be seen in Table Ql Note that even though the 
mean-field results show an almost axially symmetric so- 
lution at / = 2 (7 < 1°) the collective state shows a 
substantial 7-deformation. For the 1 = 4 state the situa- 
tion is different. We can use the pairing collective paths, 
but these are limited by the pairing collapse. Therefore 
we do not see any low lying excited states and we see no 
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FIG. 4: The wave function for the large amplitude collective 
motion in '^^Kr at I = and 2 as a function of the collective 
coordinate Q. The thick solid line is the potential energy. 
The wave functions have their scale on the right side. 



TABLE I: The deformation and excitation energy (in 
for the collective states. 


MeV) 


/n 


Oi 


O2 


O3 


2i 


22 


23 


4i 


42 


£2 

7 


0.28 
60.0 


0.37 
0.0 
0.37 


0.15 
0.0 
0.96 


0.30 
42.4 
0.32 


0.40 
9.8 
0.61 


0.36 
19.9 
1.33 


0.36 
59.3 
1.04 


0.45 
0.7 
1.07 



the prolate band-head. Our calculations support these 
conclusions with the exception of the interpretation of 
the second 0-state. We find the prolate band head to 
be lower-laying in energy then seen in experiment and 
we also see very little mixing of the oblate and prolate 
states at / = 0. An alternative interpretation would be 
that the second 0-state is a spherical state in our calcu- 
lation, since this state has a higher energy and a lager 
mixing with the oblate state. 

In summary, we have extended the method of calculat- 
ing the self-consistent collective path presented in [j, lul 
to include constraints on angular momentum. We have 
been able to determine the collective coordinate from the 
Hamiltonian without having to assume a priori which are 
the relevant degrees of freedom. The results confirm the 
importance of pairing collapsed states for the collective 
path as suggested in Ref. [111. ^^ have also seen the ef- 
fect of rotation on the collective path: Without rotation 
the path goes through a spherical saddle-point in contrast 
to the rotating case where the two minima are connected 
via the triaxial plane, and as we increase angular momen- 
tum the collective path disappears. The changes would 
probably be less pronounced if we would have allowed for 
more the one collective coordinate, which a calculation 
we intend to do in the near future. We have compared 
our calculations with experimental data |5| and found a 
reasonable agreement. 

This work was supported by the UK Engineering and 
Physical Sciences Research Council (EPSRC) 



indication of a collective wave function extended into the 
plane of non-axial deformation. 

Figure [3 shows the collective excitation spectrum for 
/ = 0, 2 and 4 in ^^Kr as a function of the angular mo- 
mentum. We can see that the prolate and oblate states 
get close in energy at finite angular momentum. We have 
also plotted the experimental data 5] for the low spin 
states. The ground state is thought to be oblate but the 
rotational band build on to of it is thought to be prolate. 
The second experimental / = state is interpreted as 
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FIG. 5: Excitation energy of the collective states for J = 0, 
2 and 4 as a function of angular momentum. Experimental 
data is taken from ^. 
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